The widths of quarkonia in quark gluon plasma 
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Recent lattice calculations showed that heavy quarkonia will survive beyond the phase transition 
temperature, and will dissolve at different temperatures depending on the type of the quarkonium. 
In this work, we calculate the thermal width of a quarkonium at finite temperature before it dissolves 
into open heavy quarks. The input of the calculation are the parton quarkonium dissociation cross 
section to NLO in QCD, the quarkonium wave function in a temperature-dependent potential from 
lattice QCD, and a thermal distribution of partons with thermal masses. We find that for the 
J/ip, the total thermal width above 1.4 T c becomes larger than 100 to 250 MeV, depending on the 
effective thermal masses of the quark and gluon, which we take between 400 to 600 MeV. Such a 
width corresponds to an effective dissociation cross section by gluons between 1.5 to 3.5 mb and 
by quarks 1 to 2 mb at 1.4 T c . However, at similar temperatures, we find a much smaller thermal 
width and effective cross section for the T. 
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INTRODUCTION 

Recently, a number of important progress have been 
made in the physics of J/tp suppression as a signature 
of quark gluon plasma that inevitably leads us to aug- 
ment the original work by Matsui and Satz[l[ with a more 
detailed study of heavy quark system at finite tempera- 
tures, before confronting the recent RHIC dataQ, and 
in predicting results for LHC. Among these theoretical 
developments are the phenomenologically successful sta- 
tistical model for J/tp production [3, 0, Hj], based on a 
coalescence assumption near T c [Bj], the recombination of 
charm pairs into J/ip 6\, and the recent lattice calcula- 
tions, showing strong evidence that the heavy quarko- 
nium will persist above T c @, 0, @, 0, EH- While these 
results seem at odd with each other, it only suggests that 
one still needs a more detailed understanding of the prop- 
erties of heavy quark system in the quark gluon plasma, 
especially between the phase transition and the dissolving 
temperatures, before a consistent picture of quarkonium 
suppression in heavy ion collision can be achieved. 

In this respect, an important quantity to investigate 
is the effective thermal width, and/or the effective dis- 
sociation cross section of a heavy quarkonium in the 
quark gluon plasma. Except for its existence, the present 
lattice results are far from making quantitative state- 
ments on the magnitude of thermal width for charmo- 
nium states above T c . Hence, in this work, we will use 
the perturbative QCD approach to calculate the ther- 
mal width. So far, such calculations have been limited 
to dissociation processes by gluons to the lowest order 



(LO) [jj, [H, [[I [H, UJ& [13], because the elementary J/ip- 
parton dissociation cross section was available only to 
that order [l^l, and to part of the next- to- leading-order 
(NLO) in the quasi free approximation ll9(. Recently, two 
of us have performed the dissociation cross section cal- 
culation to NLO in QCD[2(|. Here, we will implement 
the NLO formula, to calculate the corresponding thermal 
width of charmonium 2l| , and then perform a similar cal- 
culation for the bottonium case. 

The NLO calculation of J/^-parton dissociation cal- 
culations involves collinear divergence. When applying 
this elementary cross section to dissociation by hadrons, 
the collinear divergence is cured by mass factorization, 
which renormalizes the divergent part of the cross sec- 
tion into the parton distribution function of the hadron. 
Such complications disappear at finite temperatures, as 
the thermal masses of the partons automatically renders 
the divergence finite. The magnitude of the thermal mass 
as a function of temperature has been obtained previ- 
ously by examining the equation of state [22| . In the 
region of T c to 2T C , they are of the order of 300-400 MeV 
for quarks and 400-600 MeV for gluons. In this work, in- 
stead of following the detailed temperature dependence, 
we shall study results with a thermal mass of 400 and 
600 MeV for both the quarks and gluons. As we will see, 
with an effective thermal mass of 400 MeV, we find that 
the effective thermal dissociation cross section above 1.4 
T c is larger than 250 MeV. The NLO calculation is pro- 
portional to the derivative of the momentum space wave 
function. We will use the Coulomb wave function, whose 
size have been fitted to reproduce the result obtained by 
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FIG. 1: LO diagram 



Wong 16], using a potential extracted from lattice gauge 
thermodynamical quantities. 

In Section II, we will recapitulate the LO result. In 
Section III, we will discuss the NLO for J/ip. The LO 
and NLO results for bottonium are given in Section IV. 



T/T c 1.13 1.18 1.25 1.40 1.60 1.65 

eo(MeV) 36.4 20.9 10.1 3.4 0.14 0.004 

V< r 2 >(fm) 0.97 1.19 1.54 2.30 4.54 5.17 

oo(fm) 0.56 0.69 0.89 1.33 2.62 2.99 

TABLE I: The binding energy, the rms radius, and its corre- 
sponding Bhor radius of J/tp at finite temperature. 

perform the calculation in the rest frame of J /ip, so the 
relative velocity of J/ip and initial parton is v re i — |fe/feo|. 
The cross section is given as, 



LO RESULT 

The LO invariant matrix element for the J/tp dissoci- 
ation cross section by gluon first obtained by Peskin[l8j 
and rederived by one of us[23[ using the Bethe-Salpeter 
equation is given as, 



dp 



dpi 



— / si +7o jl -7o™ , \ 



(1) 



Here, /i, v represents the polarization index of the J/tp 
and gluon respectively, and k,pi,p2 are the four momen- 
tum of the gluon, c, c. The quantity p is the relative 
three momentum between c and c, and ip(p) is the char- 
monium wave function. M$ is the mass of a quarkonium, 
and N c is the number of color. Current conservation is 
easily shown to be satisfied, k u M^ v = 0. The energy 
conservation in the non-relativistic limit implies, 



fco + = 2m c 



|P1| 2 + \V2\ 

2m r 



(2) 



from which, the counting scheme for both the LO[l8| and 
the NLO[20| are given as follows, 



Ipll ~ |pa| ~ \p\~0(mg 2 ) 
k° ~ |fc| ~ 0(m 5 4 ). 



(3) 



The effective thermal width and cross section are ob- 
tained by folding the matrix element with the thermal 
parton distribution, n(k ) 



T eff = d 



d 3 k 



n(k Q )v re i(j(ka), 



(2tt) 3 

^-^-n(fc )er(fc )/ I 7^j"( fc o)' ( 4 ) 



(20 



where d p is the parton degeneracy, which is taken to be 
16 for the gluon in the LO calculation. Here, we will 
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where is the thermal mass of a gluon. As can be 
seen from Eq. ((5|) , the cross section is proportional to the 
absolute square of the derivative of the momentum space 
wave function, which comes from the dipole nature of the 
interaction between the gluon and the quark antiquark 
pair with opposite charges. In the calculation of the effec- 
tive thermal width or the effective cross section in Eq.(j4j, 
the cross section is integrated over the incoming energy, 
which effectively integrates over the absolute square of 
the derivative of the momentum space wave function. 
As a consequence, the results are sensitive to the size 
of the wave function only and not so much on its de- 
tailed functional form. Therefore, we will use a Coulomb 
wave function, whose Bohr radius is fitted to reproduce 
the rms radius obtained by one of us[16| by solving the 
bound states in a temperature-dependent potential ex- 
tracted from lattice gauge thermodynamical quantities. 
For the binding energy, we use the values obtained in ref. 
[l6| . Table U summarizes the J/ip binding energy, its rms 
radius, and its corresponding Bhor radius at finite tem- 
perature. The coupling constant g is set such that a s is 
0.5. As can be seen in Fig. [21 a(kn) in the LO is dominant 
around the J/ip binding energy [16j, which decreases as 
temperature increases. On the other hand, n(ko) favors 
higher fco as temperature increases. Hence the overlap in 
Eq.((4]) decreases at higher temperature, more so because 
the overlap integral starts from the thermal mass of the 
gluon. As can be seen in Fig. [3l choosing the effective 
thermal mass of the gluons to be 400 and 600 MeV, which 
is much larger than the binding energy, we find that the 
effective thermal width decreases as the temperature in- 
creases and becomes very small. With lower bound of the 
thermal mass (400 MeV) the width at LO is smaller than 
3 MeV at 1.13 T c and less than 1 MeV at 1.4 T c . At the 
NLO, there are other thermal gluon and quark induced 
interaction reactions representatively shown in Fig. 2] 
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The width due to these contributions will be calculated 
in the next section. 
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FIG. 2: (Color online) a(E g i uon ) of J/ip at LO. 
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FIG. 3: (Color online) Effective thermal width of J/ip at LO. 
The dots and circles show the results at temperatures given 
in Table U with T c = 170MeV. 



NLO RESULT 

The J/ip dissociation cross section by par tons in QCD 
at the NLO was performed by two of us[20|. The disso- 
ciation cross section can be divided into two parts; the 
dissociation by quarks and that by gluons (Fig. [4|). 

In both cases, the cross sections are given as, 



a(u) 



da 3 \M( 



(6) 



and da 3 is the 3-body phase space, 

d z k 2 d 3 px d 3 p 2 



da* = 



(27r)32fc 20 (27r)32p 10 (2irf2p 2Q 
x(27r)^(q + k 1 -k 2 - Pl -p 2 ) (7) 



Here, q, fci, k 2 , pi and p 2 are respectively the momentum 
of J/ip, incoming parton, outgoing parton, charm quark, 
and anti-charm quark. 

In the J/tp rest frame, the cross section can be written 
as. 

1 



mlm 2 ki 



dw^ 



, 2 vl — 4to 2 /u> 2 — — i 

dp A — 1 w 

16 2 7r 3 m$|fci| 



\M\ , (8) 



where p\ = (fci — k 2 ) 2 , w 2 = (q + pa) 2 - The integration 
range is 

a = (m P2 + rn pi ) 2 = Am 2 , 
13 = (u- m kl ) 2 , 
a' = -b- 



ac, 



(3' = -b+VV 2 ^ 



ac. 



(9) 



where 



b = ^ {u 2 - (m$ + m fel ) 2 } {m 2 - (m$ - ?7i fcl ) 2 } 



~2^2 ~ ( m * ~ m ki)} {™ 2 ~ m l} » 

b 2 -ac = -L{( M 2 - m 2 +m 2 J 2 - 4u 2 m^ } 

x {w 2 - (u + m kl ) 2 } {w 2 -(u- m kl ) 2 } . 

The variables appearing in a can be expressed in terms 
of w 2 ,pa and u 2 as follows, 
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ho 
k 20 
\P\ 
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V /{( W 2 -m 2 i +m|)/(2m 4 ,)} 2 
(-p 2 A + 2ml)/2, 



u 2 , 



= \ \h\ 2 



"fci' 



= fcio - (w 2 -p\ - m|)/(2ra$), 
= \J m c (k 10 - k 20 + m$ - 2m c ). 



(10) 
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Q 
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where u 2 = (q+ki) 2 , m kl is the thermal mass of a parton, 



FIG. 4: NLO diagrams induced by (a)quarks, (b)gluons. 
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FIG. 5: (Color online) a(E quark ) of q - J/ip at NLO. 



FIG. 7: (Color online) Effective thermal width of J/ip at NLO 
induced (a) by quarks and (b) by gluons 
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FIG. 6: (Color online) a(E g i uon ) of g - J/ip at NLO. 



We first consider the NLO effective thermal width in- 
duced by quark. Here, the quark degeneracy is 36 as- 
suming 3 flavors. The invariant matrix element is given 

as|20j, 



2 4 
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dp 
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Next, the gluon induced NLO calculation has the same 
degeneracy as the LO case, and the invariant matrix el- 
ement is given as follows (20j. 
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FIG. 8: (Color online) Effective cross-section of J/tp induced 
(a) by quarks (NLO) and (b) by gluons (LO+NLO) 



In the hadronic phase, the l/(fci • k 2 ) term, and the 
l/k 2 Q term give rise to collinear divergence, and soft di- 
vergence respectively. However, in QGP phase, the ther- 
mal mass of a parton in Eq. ([TTJ|) plays the role of 
a cutoff. In contrast to the LO calculation, as can be 
seen from Fig. [5] for quark and from Fig. [6] for gluon, 
a(k ) at the NLO does not vanish at large k . This is 
so because irrespective of how large the energies of the 
incoming quark or the gluon are, they can always radi- 
ate small energy gluon in the order of the binding energy 
to effectively dissociate J/ip via the LO process. Hence, 
in Eq. g]), a has non trivial overlap with the maxi- 
mum of n(kn) that increase with temperature, leading 
to the result shown in Fig. (Ha) and in Fig. [TJb) for 
quark and gluon induced NLO width respectively. One 
thing to note from Fig. [5] is that the elementary cross 
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section has a peak near threshold at 1.4 T c . Such peak 
structure only appears when the binding energy becomes 
very small and the corresponding momentum space wave 
function becomes highly peaked near zero momentum. 
When the incoming energy is small, these highly peaked 
region gives important contributions to the two dimen- 
sional phase space integral in Eq. ||5J). But when the 
incoming energy becomes large, the phase space for the 
peaked region becomes smaller and so does the total cross 
section. Such singular behavior as a function of the in- 
coming energy disappears when the binding energy be- 
comes larger. 

Here we take the thermal mass from 400 MeV to 600 
MeV[22[ within the temperature region of a few T c . With 
those masses we obtained large thermal widths. Even 
with an upper limit thermal mass of 600 MeV, the width 
exceeds 100 MeV above 1.4 T c , where we have taken 
T c =170 MeV. For example, if the thermal mass of par- 
tons is 600 MeV, and the produced J/tp remains at 1.4 
T c for 2 fm/c, its survival rate will be less than 40 %. 
As can be seen from Fig. HJa) and Fig. [8jb), with a 600 
MeV thermal mass at 1.4 T c , the effective dissociation 
cross section by a quark is about 1.0 mb and that by a 
gluon 1.5 mb. Hence, even though the J/ip might start 
forming at 1.6 T c , its effective width is very large and will 
not accumulate until the system cools down further. 



T DISSOCIATION 

Here, we present the result for the T case. The 
T wave function is less sensitive to changes in the 
temperature [3]. In the LO calculation, while the trends 
in the temperature dependence of a is similar to that 
of the J/ip, its variation in magnitudes is much smaller. 
Moreover, as can be seen in Table [H] the binding energy 
remains large. As a consequence, the overlap of a and the 
thermal distribution remains large even at large temper- 
atures, and the effective thermal width slowly decreases 
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FIG. 9: (Color online) Effective thermal width of T at LO 
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FIG. 10: (Color online) a(E quark ) of q - T at NLO. 
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FIG. 11: (Color online) a{E gluon ) of g - T at NLO. 



with temperature as can be seen in Fig. [5] 

In the NLO calculation, as shown in Fig. \TU[ and Fig. 
ITT1 the form of the cross section is similar to that of 
J/ip. However, because the binding energy of T is still 
large at high temperature, there is no peak structure near 
threshold. Moreover, because T is more tightly bound 
and has a smaller dipole size than J/ip, its corresponding 
a is also smaller. Therefore, as can be seen from Fig. 
IT27 a) for the quark induced and Fig. H2T b) for the gluon 
induced width, the overall value of the width for T is 
smaller than that of J/ip. With an upper limit of thermal 
mass of 600 MeV and at temperature of 1.65 T c , the 
sum of the LO and NLO thermal widths is less than 50 
MeV. At this temperature, the effective dissociation cross 
section by a quark is less than 0.2 mb and by a gluon 
than 0.6 mb. Therefore, unlike the J ftp case, the T has 
a smaller thermal width and effective dissociation cross 
section, and will effectively start accumulating at higher 
temperatures. Fig. H3Jb) is the sum of LO and NLO 
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SUMMARY 
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FIG. 12: (Color online) Effective thermal width of T at NLO 
induced (a) by quarks and (b) by gluons. 



In this work, we have calculated the thermal widths 
of J jt\> and T at finite temperature using the elementary 
parton-quarkonium dissociation cross section at NLO in 
QCD and assuming thermal partons with effective ther- 
mal mass. We find that for J/tp at 1.4 T c , the thermal 
width will be 100 to 250 MeV, which translates into an 
effective thermal cross section of several mb. However, 
the corrsponding width and effective cross section for the 
T is much smaller. Recently Mocsy and Petreczky [HI 
have also calculated the thermal widths of charmonium 
and bottonium in QGP, assuming that the quarkonium 
and its constituents are in thermal equilibrium with the 
surrounding. The thermal width estimated by MP is sim- 
ilar to ours for the bottonium but several times larger 
for the charmonium. The result for the charmonium by 
MP is obtained using a phenomenological formula ob- 
tained when the binding energy is much smaller than the 
temperature; hence further work has to be performed to 
understand the discrepancy. 
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FIG. 13: (Color online) Effective cross-section of T due to (a) 
a quark and (b) a gluon. 



effective cross section of T due to gluon. Because the 
contribution of LO is not negligible, the shape is quite 
different from the effective cross section due to quarks of 
Fig Ula). 



T/T c 


1.13 


1.18 


1.25 


1.4 


1.65 


1.95 


eo(MeV) 


313 


247 


203 


150 


111 


86 


V< r 2 >(fm) 


0.294 


0.331 


0.366 


0.425 


0.494 


0.562 


a (fm) 


0.17 


0.19 


0.21 


0.246 


0.285 


0.324 



TABLE II: Binding energy and rms of the T. 
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